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Abstract: We report the first Brillouin measurement of the human eye in 
vivo. We constructed a Brillouin optical scanner safe for human use by 
employing continuous-wave laser light at 780 nm at a low power of 0.7 
mW. With a single scan along the optic axis of the eye, the axial profile of 
Brillouin frequency shift was obtained with a pixel acquisition time of 0.4 s 
and axial resolution of about 60 um, showing the depth-dependent 
biomechanical properties in the cornea and lens. 
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1. Introduction 

Brillouin optical microscopy is a novel technique based on Brillouin light scattering 
spectroscopy for probing the viscoelastic properties of a sample with three-dimensional 
resolution [1]. Brillouin spectroscopy has long been applied for material characterization and 
environmental sensing [2, 3]. The first Brillouin spectroscopy of biological tissues was 
demonstrated in the early 1980s [4, 5], where the Brillouin spectra of lens and cornea at single 
spatial points were obtained with an integration time of 10 min to 1 hour. The recent 
development of parallel spectrometer drastically reduced the acquisition time to 1 s or less [6], 
which allowed spatially resolved measurements [1,7] and imaging [8]. An in vivo study of 
murine lenses revealed age-related changes in the Brillouin profiles of the lens nucleus and 
also found an empirical relationship between the Brillouin modulus and conventional Young's 
modulus [8]. Spatially-resolved measurements of human lenses [9] and animal corneas ex vivo 
[10] were recently demonstrated. 

The possibility of applying the technique to a living patient offers an exciting opportunity 
for early detection of ocular problems, such as presbyopia and corneal ectasia, and for 
understanding their underlying mechanisms [11]. Here, we report the development of a 
clinically viable Brillouin microscope and describe the first Brillouin measurement of the 
human eye in vivo. For human use, the instrument employs a low-power laser light at 780 nm 
and a parallel Brillouin spectrometer optimized for the infrared wavelength. 

2. Setup and methods 

Figure 1 shows a schematic of the experimental setup. The instrument is comprised of three 
parts: laser- scanning confocal microscope, human interface, and VIPA-etalon spectrometer. 
The light source was a continuous-wave (cw) external-cavity semiconductor laser emitting a 
single longitudinal mode at 780 nm (Toptica Photonics). The laser output was attenuated with 
a neutral density filter and expanded to the 1/e 2 laser beam diameter of about 7 mm. A 10/90 
beam sampler directed 10% of the laser power to the human interface through a beam- 
scanning unit. We used two mirrors for x-y transverse beam positioning and a motorized 
translation stage (max speed 1.5 cm/s, Zaber Technologies) for z-scan. We used two different 
objective lenses: an achromat (Thorlabs, f = 35 mm, NA = 0.055) with a long working 
distance (35 mm) was used; and a 4X microscope objective (Edmund Optics, NA = 0.1) with 
higher spatial resolution was utilized. The backward scattered light from the eye was collected 
with a single-mode optical fiber through the beam sampler (90% transmission). We used two- 
stage VIPA etalons optimized for the near-infrared wavelength (Rl = 99.9%, R2 = 98% at 
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760-800 nm), configured with the cross-axis cascade principle [12]. In each stage, the light 
beam was line-focused into the input window of the VIPA etalon tilted by about 2 degrees 
and the output pattern went through a spatial filter to remove the elastic scattering light 
component. The spectrally dispersed pattern was imaged onto an EM-CCD camera (Andor, 
IXon DU-897). The camera was operated at full gain with a negligible readout noise. 




Fig. 1. Schematic of the Brillouin confocal microscope. OI: Optical isolator. ND: Neutral 
density filter wheel; BS: beam sampler; SMF: single mode fiber. 

The human interface was a modified version of standard ophthalmic slit-lamp instrument. 
It has a chin support and a headrest, which allowed the height of the head of the human 
subject to be adjusted with a manual translation stage. All measurements were performed in 
the relaxed accommodation state. Prior to the Brillouin scan, we aligned the position of the 
beam to the corneal apex by finding the location that maximizes the optical back-reflection 
from the corneal surface. The focus of the laser beam was scanned along the optical axis from 
the cornea to the lens or backward. Axial scanning was performed by moving the objective 
lens with speeds ranging between 30 and 80 um/s. The Brillouin spectra were recorded in the 
CCD camera with the frame integration times between 0.2 and 1 s. The laser power at the 
sample could be adjusted between 0.7 and 3mW. Beam scanning, CCD acquisition and 
spectral analysis were synchronized with LAB VIEW and MATLAB programs. The 
measurement results were displayed in real-time on a computer. The axial profile of the 
Brillouin shift was plotted with respect to the focal depth in the eye, calculated using the 
average refractive index (n = 1 .4) of the eye anterior chamber, which was measured by a low- 
coherence interferometer. 

3. Safety consideration 

We considered radiation-induced damage thresholds in the retina, as well the cornea and lens. 
Maximum permissible exposure (MPE) is defined as the highest power or energy density that 
can be admitted to the eye without causing a biohazard. MPE corresponds to 10% of the dose 
that has a 50% chance of creating damage in worst-case scenario conditions. In the Brillouin 
scanner, the laser light is focused in the anterior segment and diffused onto the retina (Fig. 2). 
According to the International Commission on Non-Ionizing Radiation Protection (ICNIRP), 
for cw sources in the wavelength region 400-1050 nm the exposure limit for cornea-lens 
thermal safety is 4W/cm 2 ; i.e. MPE = 32 mW in a 1-mm-diameter zone (0.79 mm 2 in area) 
[13]. The 1-mm aperture for irradiance averaging is based on the thermal modeling [14], 
which shows that the temperature rise is independent of the beam size up to 1 mm due to 
rapid thermal conduction (-0.5 W/m/°C) in the cornea and lens. As for the retina, the 
dominant mechanism of retinal damage is thermal for an exposure time (T) longer than 0.25 s. 
Considering the distance 17 mm between the lens and retina and a = 0.2, i.e. NA = 0.1, the 
beam size on the retina is >3 mm. Applying the same thermal limit of 4 W/cm 2 ,we calculate 
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MPE = -300 mW. According to the guidelines from American National Standard Institute, 
the exposure limit for the retinal thermal hazard has been expressed as 1.8xl0~ 3 C A C E T 025 
[W/cm 2 ], where C A = 1.45 (for = 780 nm), C E = 267 (for a = 0.2; i.e. NA = 0.1), and the 
aperture size of 7mm (area of 0.38 cm 2 ) [15]. The MPE is calculated to be 265 mW for T = 1 
sand 95 mW for T = 60 s, which is consistent with the above analysis. Taken together, we 
conclude that the Brillouin scan in the cornea and lens does not pose risk to the human eye. 




Fig. 2. Laser illumination in the Brillouin scanner. By focusing the beam in the transparent 
tissue of the anterior chamber (cornea, aqueous humor and lens), the beam is diverging and 
covers a large area when it hits the retina. 

4. Results 

The experimental protocol has been reviewed and approved by the Institutional Review Board 
for human research at the Massachusetts General Hospital. A 42-year-old male healthy 
volunteer was enrolled. Prior to the Brillouin test, a full eye exam including topography and 
pachymetry was conducted, which did not reveal any abnormality and measured a central 
corneal thickness of 530 um. The optical power used in the Brillouin scan was measured to be 
0.7 mW at the eye, about 50 times lower than the MPE. 




Brillouin shift (GHz) 

Fig. 3. Representative Brillouin spectra (anti-Stokes peaks) from a human eye. (a) Corneal 
stroma, (b) Aqueous humor, (c) Lens nucleus, (d) Vitreous humor. Top panels: raw CCD 
spectra (average of four frames). Bottom panels: corresponding spectra. 

Figure 3 shows representative spectra acquired in four regions of the eye: cornea, aqueous 
humor, lens nucleus and vitreous humor, respectively. The top panels are raw CCD images 
featuring anti-Stokes Brillouin scattering light, and the bottom panels are the corresponding 
Brillouin spectra extracted from the CCD frames by the image analysis that has been 
described elsewhere [6,8]. Briefly, a reference arm (not shown in Fig. 1) periodically 
interrogates samples, such as water and PMMA, of known Brillouin frequency shifts. From 
the reference spectra, the direction of the spectral dispersion axis (diagonal lines) and its 
pixel-to-GHz ratio are measured (0.38 GHz/pixel). To obtain the spectrum along the spectral 
dispersion axis, the CCD images are interpolated (x4) and integrated over the orthogonal axis. 
The Brillouin frequency shift is determined with Lorentzian curve fitting. The measurement 
sensitivity of the Brillouin peak detection was measured to be 30-120 MHz for this data set. 
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Figure 4 shows the measured axial profile obtained from the left eye of the volunteer. For 
Brillouin scan, the beam focus was translated along the optic axis of the eye at a speed of 50 
um/s (in air), and the Brillouin spectra were acquired with a CCD frame rate of 2.5 Hz (i.e. 
integration time of 0.4 s). The axial profile of the lens shown in Fig. 4(a) was obtained with 
the 35-mm achromatic lens with an axial resolution of about 350 um. It features a typical bell 
shape with a slope in the anterior cortex, a plateau at about 6.05 GHz in the lens nucleus, and 
a decline in the posterior cortex towards the vitreous humor. The axial profile shown in 
Fig. 4(b) was obtained with the 4X objective at a higher axial resolution of -60 um. It reveals 
a better resolved anterior cortex of the lens, with the Brillouin shift varying steeply from 5.25 
GHz to 5.85 GHz in the initial 300 um layer. The high-resolution scan also allowed us to 
obtain the Brillouin profile of the cornea. Between 100 um and 400 um depths, corresponding 
to the corneal stroma, the Brillouin shift declined slowly from -5.6 GHz to -5.5 GHz. In the 
posterior region of the cornea, the Brillouin profile showed a steeper decline from 5.5 GHz to 
5.25 GHz in about 200 um. With the finite confocal resolution, strong Fresnel reflection from 
the corneal surface made it difficult to analyze the first 70 um from the corneal surface. The 
aqueous humor was measured to have a fairly constant Brillouin frequency, slightly higher 
than the Brillouin shift of pure water. 
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Fig. 4. In vivo Brillouin measurement of a human eye. (a) Depth profile of the crystalline lens, 
(b) Depth profile of the anterior chamber. The data are an average of four such scans taken ten 
minutes apart. 

We performed a total of four scans on the same subject within the span of one month. The 
variation of the Brillouin shift in the aqueous humor was measured to be about 30 MHz. As 
per approved protocol, at eight months after the Brillouin tests the subject underwent a 
follow-up full eye exam. As expected, the results showed no abnormalities and no changes in 
static or dynamic visual acuity. 

6. Discussion 

The measured Brillouin frequency shifts of the human lens and cornea in vivo agreed well 
with previous measurements of ex vivo tissues. Vaughan and Randall reported Brillouin shifts 
of 4.54 GHz, 4.71 GHz, and 4.71 GHz for cornea, cortex and nucleus, respectively [4, 5]. 
Bailey et al. reported Brillouin shifts of 7.75 GHz (cortex) and 8.8 GHz (nucleus) [9]. 

Brillouin shift, £2, can be written as Q = 2n/ A^M I p cos(# / 2) where n, M, and p are the 

refractive index, longitudinal elastic modulus (hereinafter "Brillouin" modulus), and density, 
respectively, X is the wavelength of probe light, and 6 is the angle between incident and 
scattered light. For comparison, we converted the Brillouin shifts to Brillouin modulus M 
using n = 1.37, 1.37 and 1.42, and p = 1.11, 1.06 and 1.11 g/ml for cornea, cortex and nucleus, 
respectively [4, 16]. For cornea, the modulus is computed to be 2.7 GPa [4]; for cortex, the 
moduli are 2.8 GPa [4, 5], and 2.36 + 0.09 GPa [9]; for nucleus, 2.93 GPa [4, 5], and 2.79 + 
0.14 GPa [9]. In our measurements corneal Brillouin modulus ranged from 2.82 to 2.5 GPa 
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and the lens Brillouin moduli varied between 2.38 and 3.1 GPa. These values are consistent 
with the values reported by Vaughan and Bailey et al [5, 9]. One caveat in this analysis is that 
literature data on index and density are variable, in part because index and density vary within 
the ocular tissues [16, 17]. Thus, the conversion to Brillouin modulus has an uncertainty by 2- 
3%. 

Our data showed a marked depth dependence of longitudinal modulus in the cornea and 
lens. Previous mechanical tests hinted at elasticity gradient within both cornea and lens tissue. 
Kohlhaas et al. measured the anterior corneal stroma to have three times higher Young's 
modulus than the posterior stroma [18]. Randleman et al. reported that the elastic modulus 
decreased from the front to the back of the cornea [19], consistent with our in vivo data and 
previous animal data [10]. As for the lens, the crystalline proteins are packed with 
increasingly tighter density toward the nucleus so that a varying modulus is expected. The 
increase in density of crystalline proteins was measured in animals [20] . The spatial variation 
of Brillouin modulus has been reported in the lens in animals ex vivo [7] and in vivo [8]. 

Brillouin modulus differs from quasi-static Young's modulus from a fundamental point of 
view. For a biological tissue with high water content, the Brillouin modulus is much higher 
than Young's modulus. This effect produced an apparently counterintuitive result: while the 
corneal tissue is known to have higher Young's modulus than the lens tissue, the Brillouin 
modulus of the cornea is actually lower mostly because of the higher water content in the 
cornea (-80%) than in the lens (-65%). Nevertheless, a good correlation can be found 
between Brillouin modulus and Young's modulus for the same tissue type [8]. 

The biomechanical properties of the lens and cornea have been linked to various ocular 
disorders. Presbyopia is thought to occur due to the age-related stiffening of the lens. 
Although Brillouin microscopy was shown to be able to detect age-related stiffening in 
animals [8], its ability to detect this effect in human lenses has been questioned [9]. Stiffening 
of the human lens has also been linked to the onset of age-related nuclear cataracts [21]. As 
for the cornea, keratoconus and corneal ectasia are associated to a reduced corneal modulus 
[22] . However, current clinical analysis of the eye is limited to structural measurements while 
traditional mechanical tests are destructive and cannot be employed in clinical settings. The 
only instrument for corneal mechanical analysis is the ocular response analyzer (ORA) [23], 
but its sensitivity and clinical usefulness are still under question [24] . 

The spatial resolution of the current Brillouin scanner does not allow us to single out the 
Brillouin scattering signature of the lens capsule. The contribution of the capsule to the 
accommodative amplitude has been studied, although the mechanical modulus of the capsule 
may not change much with age [25]. Using a higher-NA objective lens, it may be possible to 
resolve the capsule and measure its Brillouin modulus. Another limitation of current Brillouin 
instrument for use in the clinic is the relatively long acquisition time. A single full axial scan 
across the eye takes about 1 min, although the scan time could be reduced to <10 s at the 
expense of the scan range or sampling intervals. By reducing the optical loss in the 
spectrometer and improving VIPA etalon coating procedures, it might be possible to decrease 
the spectral acquisition time by 10 fold and thereby to obtain several axial or transverse scans 
across the lens. This speed improvement will help to assess the equatorial distribution of the 
elastic modulus of the lens. 
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